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' Abstract 

^ \ We show that globally and regularly hyperbolic future geodesically incomplete 



X 



isotropic universes, except for the standard all-encompassing 'big crunch', can ac- 
commodate singularities of only one kind, namely, those having a non-integrable 



^ \ Hubble parameter, H. We analyze several examples from recent literature which 

illustrate this result and show that such behaviour may arise in a number of dif- 
ferent ways. We also discuss the existence of new types of lapse singularities in 
inhomogeneous models, impossible to meet in the isotropic ones. 



1 Introduction 



It has been pointed out very recently that finite-time, finite Hubble parameter singular- 
ities can occur in the expanding direction of even the simplest FRW universes in general 
relativity and other metric theories of gravity and, although such behaviour depends on 
the details of the particular model, yet, it is met in a very wide and extremely varied 
collection of possible cosmologies (cf. Refs P-jSO] for a partial list). These singulari- 
ties cannot be accommodated by the usual singularity theorems for their features 
arise from necessary, not sufficient, conditions for their existence while their character is 
somewhat milder than the standard, all-encompassing 'big crunch' singularities met in 
other model universes. This situation appears somewhat wanting as all known relevant 
models provide merely examples, not general theorems, of such a phenomenon. 

What is the underlying reason for such a behaviour? In this paper, we show that the 
finite-time, finite- if singularities met in isotropic cosmologies satisfy all the assumptions 
of the completeness theorem of pSj except the one about an infinite proper-time interval 
of existence of privileged observers and because of this reason are geodesically incomplete. 
This leads to a classification of the possible types of singularities in the isotropic category 
based on necessary analytic conditions for their existence. We also briefiy discuss how 
new types of such singularities may arise in more general inhomogeneous models. 

2 Necessary conditions for future singularity forma- 
tion 

Although we shall focus exclusively on isotropic models, it is instructive to begin our 
analysis by taking a more general stance. Consider a spacetime (V, g) with V = Ai x I, 
I = {tQ,oo), where is a smooth manifold of dimension n and ("+^)(^ a Lorentzian 
metric which in the usual n + 1 splitting, reads 
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Here N = N{t,x'^) is called the lapse function, (3^{t,x^) is called the shift function and 
the spatial slices Ait {= A4 y< {t}) are spacelike submanifolds endowed with the time- 
dependent spatial metric gt = gijdx^dxK We call such a spacetime a sliced space jHSI- A 
sliced space is time-oriented by increasing t and we choose I = {to, oo) because we shall 
study the future singularity behaviour of an expanding universe with a singularity in the 
past, for instance at t = < to- However, since t is just a coordinate, our study could 
apply as well to any interval X C M. 

A natural causal assumption for (V, g) is that it is globally hyperbolic. This implies 
the existence of a time function on {V,g). In a globally hyperbolic space the set of all 
timelike paths joining two points is compact in the set of paths, and spacetime splits as 
above with each spacelike slice Ait a Cauchy surface, i.e., such that each timelike and 
null path without end points cuts Ait exactly once j31j. 

We say that a sliced space has uniformly bounded lapse if the lapse function N is 
bounded below and above by positive numbers and Nm, 

0<Nm<N< Nm. (2.2) 

A sliced space has uniformly bounded shift if the gt norm of the shift vector (3, projection 
on the tangent space to TVIj of the tangent to the lines {x} x X, is uniformly bounded 
by a number B. 

A sliced space has uniformly bounded spatial metric if the time-dependent metric 
gt = gijdx^dx^ is uniformly bounded below for alH G X by a metric 'J = gto, that is there 
exists a number A > such that for all tangent vectors f to it holds that 

A^.jvV < g.jv'v^. (2.3) 

A sliced space (V, g) with uniformly bounded lapse, shift and spatial metric is called 
regularly hyperbolic. 

Denoting by VA^ the space gradient of the lapse A^, by Kij = — A^r°- the extrinsic 
curvature of A^f, and by \K\^ the product g°-^g^^KabKij, we have the following theorem 
which gives sufficient conditions for geodesic completeness jH^ : 
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Theorem 2.1 Let {V,g) be a sliced space such that the following assumptions hold: 
CI (y,g) is globally hyperbolic 
C2 (V, g) is regularly hyperbolic 

C3 For each finite ti, the space gradient of the lapse, |VA^|g, is bounded by a function 
oft which is integrable on [ti,+oo) 

C4 For each finite ti, \K\g is bounded by a function oft which is integrable on [ti, +00). 

Then (V, g) is future timelike and null geodesically complete. 

It is known that in a regularly hyperbolic spacetime, condition CI is in fact equivalent 
to the condition that each slice of {V,g) is a complete Riemannian manifold, cf. [33J. 
The completeness theorem gives sufficient conditions for timelike and null geodesic 
completeness and therefore implies that the negations of each one of conditions CI — CA 
are necessary conditions for the existence of singularities (while, the singularity theorems 
(cf. [31 ) provide sufficient conditions for this purpose). Thus Theorem 12.11 is precisely 
what is needed for the analysis of models with a big rip singularity inasmuch these models 
require necessary conditions for future singularities. 

Now for an FRW metric, = 1 and /? = 0, and so condition C2 is satisfied provided 
the scale factor a{t) is a bounded from below function of the proper time t on X. Condition 
C3 is trivially satisfied and this is again true in more general homogeneous cosmologies 
where the lapse function N depends only on the time and not on the space variables. 
In all these cases its space gradient, |VA^|g, is zero. Condition C4 is the only other 
condition which can create a problem. For an isotropic universe \K\g' = 3(d/a)2 = 3H^, 
and so we conclude that FRW universes in which the scale factor is bounded below can 
fail to be complete only when C4 is violated. This can happen in only one way: There is 
a finite time ti for which H fails to be integrable on the time interval [ti, 00). Since this 
non-integrability of H can be implemented in different ways, we arrive at the following 
result for the types of future singularities that can occur in isotropic universes. 
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Theorem 2.2 Necessary conditions for the existence of future singularities in globally 
hyperbolic, regularly hyperbolic FRW universes are: 

51 For each finite t, H is non-integrable on [ti,t], or 

52 H blows up in a finite time, or 

53 H is defined and integrable (that is bounded, finite) for only a finite proper time 

interval. 

When does Condition 5*1 hold? It is weU known that a function H{t) is integrable 
on an interval [ti,t] if H{t) is defined on [ti,t], is continuous on (ti,t) and the limits 
lim^_^^+ H{t) and lim^^j- H{t) exist. Therefore there are a number of different ways 
which can lead to a singularity of the type SI and such singularities are in a sense more 
subtle than the usual ones predicted by the singularity theorems. For instance, they may 
correspond to 'sudden' singularities (see [19j for this terminology) located at the right 
end (say tg) at which H is defined and finite but the left limit, lim^^^+ H{t), may fail 
to exist, thus making H non-integrable on [ti,ts], for any finite ts (which is of course 
arbitrary but fixed from the start). We shall see examples of this behaviour in the next 
Section. Condition 5*2 leads to what is called here a blow-up singularity corresponding to 
a future singularity characterized by a blow-up in the Hubble parameter. Note that SI 
is not implied by 5*2 for if H blows up at some finite time tg after ti, then it may still be 
integrable on [ti,t], ti < t < tg. Condition 5*3 also may lead to a big-rip type singularity, 
but for this to be a genuine type of singularity (in the sense of geodesic incompleteness) 
one needs to demonstrate that the metric is non-extendible to a larger interval. 



3 Examples 

We give below some representative examples to illustrate the results of the previous 
Section. An example of a blow-up singularity is the recollapsing flat FRW model filled 
with dust and a scalar field with a 'multiple' exponential potential considered in jTHj . 
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These authors choose the potential to be of the form V{(j)) = Wq—Vq sinh(y^3/2/t0) where 
Wq and Vo are arbitrary constants and k, = ^/SttG/^, and spht the scale factor according 
to the transformation = xy. Here ones sets x = C[expxicosx2 + exp(— xi) C0SX2], 
y = C[expxisinx2 + exp(-xi)sin(-X2)] with xi = wi{t-to),X2 = W2{t-To), where 
C > 0, to an arbitrary constant, Tq is the 'initial' time and Wi, W2 positive parameters 
such that wi — W2 = 3/4fi;^PVo! 2^11^2 = 3/4fi;^Vo. Using the results of Section 2, we 
can immediately see why this model has such a singularity (in ^H] this was proved 
by different methods). For large and positive values of the time parameter t the scale 
factor becomes a = C^/^ exp(2/3xi) cos^/^ X2 sin^^'^ X2, which is obviously divergent. We 
then find \K\g^ to be essentially (that is except unimportant constants) proportional to 
W2(cotx2 — tanx2)/3, and so the extrinsic curvature blows up as t tends to the finite 
time value Tr/{2w2) + Tq. This is a blow-up singularity in the sense of Condition 52 of 
Section 2. 

Another example of a future blow-up singularity is provided by the 'phantom' cos- 
mologies (see relevant references in [^). In all these models different methods are used, 
depending each time on the analysis of the field equations of each phantom model, to 
prove the existence of such future singularities. For instance, in jTTj we meet a flat FRW 
universe with a minimally coupled scalar field (p in Einstein's gravity with the equation 
of state p = wp, with w < —1 (phantom dark energy). In this model the scale factor 
takes the form 



ao^(^-^-)/^ + ^^i±|^(t-to) 



3{l+tu) 

(3.1) 



where A = 8ttGC/3 and C an integration constant. Hence the extrinsic curvature 
becomes 



\K\l = 3A 



(3.2) 



and so when w < —1 and t = to + 2/[3y/A{\w\ — l)]^^^ ^^le extrinsic curvature 

diverges thus causing a blow-up singularity. Similar behaviour is found in other singular 
phantom cosmologies, for instance those of Refs. [I], 0, [12]; IIS]; [21]; namely, there is 
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some finite time t/ at which the Hubble parameter blows up. 

Had the evolution was characterized by an integrable H for every t, then we would 
expect from Theorem 12. II all these models to be timelike and null geodesically complete. 
An example illustrating this is given by the flat FRW universe filled with phantom dark 
energy which behaves simultaneously as Chaplygin gas studied in |25] (see also |I3j). 
Here the phantom dark energy component satisfies the equation of state p = wp with 
w < —1, as well as the equation of state of a Chaplygin gas p = —A/p, where A is a 
positive constant (similar results will hold when w G (—1,-1/3) -fc-essence models, see, 
e.g., [Q). Then integrating the continuity equation p = — 3a(p + p)/a, we find 

p\t)=A + {pl-A){ao/a{t)f. (3.3) 

or, using the two equations of state and substituting to ()3.3|) leads to 

p(t) = po'h^^o + (1 + Wo){ao/a{t)ff' (3.4) 

where A = —wopo^, Wq < —1. Finally, the solution for the scale factor reads a{t) = 
(^gCi(i-io)^ and so the Hubble parameter in this case is constant, H = Ci/6. This non 
singular behavior is due to the fact that all the conditions of the Completeness Theorem 
12. II are met. 

We now move on to an example of a big rip singularity in the sense Condition 5*1 
of Section 2. Barrow in ^H] rightly calls such singularities 'sudden'. He considers an 
expanding FRW model with a fluid that satisfies the energy conditions p > and p+3p > 
0, and shows that the special solution 

a(t) = 1 + (^^y (a. - 1) + (1 - ^)^ (3.5) 

with l<n<2,0<g<l and = a{ts), exists, as a smooth solution, only on the 
interval {0,ts), while and Hg = H{ts) are finite at the right end. Note that setting 
A = q{as - l)/t1,B = n/t"^, we find that d = At"-'^ + B{ts - t)"~^ which means that, 
unless q = 1, a blows up as t — 0, making H continuous only at (0, t^). Also a(0) is finite 
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and we can extend H and define it to be finite also at 0, H{Q) = Hq, so that H is defined 
on [0,ts]). However, since \imt^Q+ H{t) = ±00, we conclude that this model universe 
implements exactly Condition 5*1 of the previous Section and thus H is non-integrable 
on [0,ts], ts arbitrary. This is a big rip singularity characterized by the fact that as 
t tg one obtains a —00. Then using the field equation we see that this is really a 
divergence in the pressure, p/2 —a/a — H'^/2 — k/{2a^), and so p — > 00. In particular, 
we cannot have in this universe a family of privileged observers each having an infinite 
proper time and finite H. 

Another way to see this result in the particular case when p > and p > 0, is to use 
the following result from pp. 141-2: In an FRW universe with p > and p > 0, 
given any vector X at any point g, the geodesic 7(1^) which passes through the point 
g = 7(0) in the direction of X, is such that either 

• 7(1^) can be extended to arbitrary values of its affine parameter or 

• there is some value I'o > such that the scalar invariant {Rap — ^/'^Rga/3){R°''^ — 
l/2Rg"f^) is unbounded on 7([0,i;]). 

Therefore under the assumptions of [11], the invariant {Rap — l/2Rgai3){R°'^ — l/2Rg°'^) 
is calculated to be 

12 24a2 12a^ 12a Ua^d VId? , , 

— + ^ + ^ + — + —T + (3.6) 

or or or a-^ a'^ 

and since a a{ts), H{t) Hs, p{t) 00, d/a —00 as t ^ ts, we see that 
{Ra/3 — ^/'^Rgai3){R°'^ + l/2Rg°''^) is unbounded at tg- Hence we find that this spacetime 
is geodesically incomplete. 

Further, Carroll et al in [H] provide a different model leading, however, to a similar 
future big rip singularity. They start from a theory of gravity with lagrangian density 
L = R+aR~'^ + Lmat and conformally transform to an Einstein frame to get an equivalent 
theory which is described as general relativity coupled to a system comprised by the 
scalar field resulting from the conformal transformation and the conformally transformed 
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matter. Since the two matter components in the Einstein frame are not separately 
conserved, there is now a non-trivial scalar field-matter interaction which manifests itself 
with new terms appearing in the field equations. It then follows that as — 0, the time 
derivative of the extrinsic curvature blows up, H oo, while H itself is finite at = 
(cf. [14]). This is clearly a big-rip type singularity. 

A different example of a big rip singularity is given by Borde et aim I3j . They showed 
that in infiating spacetimes just a bounded averaged-out Hubble function is necessary to 
produce a singularity with finite H . Why do they obtain such a behaviour? To answer 
this question we again use our Theorem 12.21 and conclude that their condition given 
by Eq. (11) in ^ is precisely the negation of condition C4 of Theorem 12.11 above. In 
effect, what these authors do is to find that a necessary condition for a singularity in 
the spacetimes considered is that H be finite on a finite time interval. One recognizes 
that under their assumptions, privileged observers cannot exist in infiating universes for 
an infinite proper time, because had such observers existed. Theorem 12.11 would then 
imply that these infiating spacetimes are timelike and null geodesically complete. More 
precisely, in |3j one starts with a fiat FRW model, ds^ = dt^ — o?{t)dx^ , and considers all 
quantities along a null geodesic with affine parameter A. Since the model is conformally 
Minkowski we have d\ oc a(t)dt, or, dX = a{t)dt/a{ts), so that dX/dt = 1 for t = ts, 
where ts is a finite value of time. One then defines 

Ha. = . / H{X)dX, (3.7) 

and so if Hav > one finds 

~ X{ts) - A(tO A(,^) ~ A(t.) - A(t.) ^ - X{ts) - X{t,y 

(3.8) 

This shows that the affine parameter must take values only in a finite interval which 
implies geodesic incompleteness. A similar proof is obtained for the case of a timelike 
geodesic. Notice that condition ()3.8j) holds if and only if condition C4 of Theorem 
12. II is valid only for a finite interval of time, thus leading to incompleteness according to 



8 



Theorem 12 .21 A similar bound for the Hubble parameter is obtained in j^j for the general 
case, and one therefore concludes that such a model must be geodesically incomplete. 

A relaxation of the requirement that H be finite for only a finite amount of proper 
time leads, as expected by Theorem 12 .1^ to singularity-free inflationary models evading 
the previously encountered singularity behaviour. Such models have been considered 
in J3] (see also |H])- It is particularly interesting in our context to see why such 
behaviour actually occurs: In [E] , the scale factor is assumed to be bounded below by a 
constant (so that this universe is regularly hyperbolic in our notation), and its form is 
given by (these authors also take this to be qualitatively true in more general models). 



One then finds that H is not only bounded as t — > — oo, but actually becomes zero 
asymptotically. Thus by Theorem 12.11 this universe is geodesically complete. It is in- 
structive, however, to notice that a slight change in the inflationary solution can drasti- 
cally alter this behaviour. Consider for example the old quasi-exponential models first 
considered by Starobinski (see and references therein), where the scale factor is given 
by a = Co exp [cit + C2t^], with the q's constants. Then we see that H = ci + 2c2t and so 
it blows up to ±oo as t ^ ±oo depending on the sign of the constant C2- It is interesting 
to further notice (cf. jHH]) that this solution is an attractor of all homogeneous, isotropic 
solutions of higher order gravity theories and so we may conclude that the blow up be- 
haviour leading to a singular regime is perhaps a more generic feature than completeness 
in such contexts. 

Examples of big rip as well as geodesically complete types of evolution are also met in 
many brane cosmologies, cf. |E], [Zj, [IDj, [H], [2E]- The singularities in such models 
are all characterized by the fact that there exists no admissible slicing with an infinite 
proper time of privileged observers and the Hubble parameter remains finite but only 
for this finite interval of proper time and cannot be defined beyond that interval (cf. 
[Sj), and hence in accordance with Theorem 12.21 these spacetimes must be geodesically 
incomplete. 




(3.9) 
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However, geodesically complete solutions also exist and such phases in brane evolution 
are described in detail by the series of brane models studied in The main idea is 
rather simple, cf. [7j, [22] • We start with a fiat 3-brane filled with ordinary matter as well 
as phantom dark energy embedded in a 5- dimensional bulk. This model goes through 
a series of cycles of finite accelerating expansion and contraction. At the end of each 
cycle it bounces avoiding in this way a singular behavior. In the expansion phase the 
phantom dark energy component increases and drives cosmic acceleration. This rapid 
acceleration tears apart any bound structure produced during expansion and towards the 
end of the expanding phase phantom dark energy has become so high that modifications 
of the Friedman equation become important. The modified Friedman equation on the 
3-brane is 

A4 Svr , 47r ,9 9 C ^ 
if =— H ^P + e( ^) p + -r, (3.10) 

where C is an integration constant, e corresponds to the metric signature of the extra 
dimension and A4 is the cosmological constant on the brane. Assuming that C/a^ is 
negligible and e < 0, we have for the critical brane case (A4=0) 

f^2^^(p_4L), (3.11) 

hence, H = when pbounce = 2|cr| where cr is the tension of the brane. At this scale 
the model will turn around and start to contract. During contraction the phantom 
component will decrease but at late times, when the model becomes matter or radiation 
dominated, the energy density will be very high and modifications of the Friedman 
equation will again become important. When the critical value of the density is reached 
the model will bounce and start to expand. Using the assumptions of this model, we 
see from Eq. ()3.1H1 that the Hubble parameter remains eternally finite and also all 
other assumptions of Theorem 12.11 are satisfied. Therefore this universe is geodesically 
complete. 
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4 Discussion 



The main result of this paper is the recognition that future big rip singularities occur 
because there is a finite time, say at 0, such that the Hubble parameter is not integrable 
on [0, oo). In such universes privileged observers cannot exist for an infinite proper time 
starting from 0. We know from Theorem 12.11 (cf. that if such observers existed 

for an infinite proper time, then such a universe would be timelike and null geodesically 
complete. This condition is not satisfied in recent models with big rip singularities and 
therefore such universes may not be complete. In fact, examples show, cf. Refs. jS], |14j . 
[T^ . that such spaces are necessarily singular. 

It is interesting that isotropic models coming from completely different motivations 
reveal their tendencies to have future singularities (that is future geodesic incomplete- 
ness) for exactly one the same reason, namely the non-integrability of H. Of course, the 
fact that this condition can be sustained in different models stems from the different ways 
these models are constructed and the specific features they share. It may be supported 
by the presence of a scalar field inducing a novel interaction with matter as in the case 
of f{R) theory considered in Ref. Jl], or lead to a divergence in the fluid pressure and a 
corresponding one in a as in jTHj, or in the inflationary character of the particular model 
as in 3^, etc. 

We have further shown that when we restrict attention to homogeneous and isotropic 
universes, these non-integrable-if singularities are the only types of future big-rip sin- 
gularities which occur. This leads us to the following question: Is it possible to find a 
globally hyperbolic, regularly hyperbolic inhomogeneous cosmology which satisfies con- 
dition C4 of Theorem 12.11 but not C3? Since |VA^|g is now a function not only of the 
time but also of the space variables it is not automatically zero. If true, this effect will 
lead, in addition to the ones already present in the isotropic case (i.e., those having 
an integrable |VA^|g but failing C4), to two new types of lapse singularities both not 
satisfying condition C3. Namely, that which has a diverging |VA^|g in a finite time cor- 
responding to a blow-up lapse singularity and secondly, that with |VA^|c, finite only for 
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a finite interval of proper time, a big-rip lapse singularity. The latter singularity, if it 
exists, will necessarily have several qualitative features distinct from the corresponding 
blow-up or big-rip (extrinsic) curvature singularities discussed here. In turn, these will 
not anymore be the only available types of singularity in the inhomogeneous category, 
as indeed they are for isotropic universes. 
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